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ABSTRACT 


Theoretical  single -station  and  network  signal  detection  capability 

models  for  surface  waves  are  given  in  this  report.  Application  of  these  models 

to  world-wide  signal  detection  capability  by  the  Very  Long  Period  Experiment 

stations  and  the  ALP  A,  NORSAR,  LASA  arrays  show  high  probability  of  at  least 

2 -station  detection  at  about  M  =  3.  0  _  3. 1  for  most  of  Eurasia  for  both  20  and 

8 

40- second  Rayleigh  waves.  More  rapid  attenuation  of  the  40  second  wave  ampli¬ 
tudes  with  distance  compared  to  20  second  amplitudes  is  demonstrated  and  in¬ 
cluded  in  the  theoretical  capability  model. 

Evaluation  of  a  network  of  such  stations  to  decrease  problems  in 
detection,  unmixing  of  signals,  and  evaluating  source  radiation  pattern  effects 
is  discussed.  Preliminary  data  for  beamsteering  continental  and  world  wide 
long  period  array  networks  is  inconclusive  as  yet  because  of  a  limited  multiple 
station  observation  base. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  necessarily  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical 
Applications  Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 

This  report  describes  investigations  of  the  network  aspects  of  the 
Very  Long  Period  Experiment  (VLPE)  stations,  with  the  purpose  of  developing  a 
basis  for  identification  of  both  the  strong  and  potentially  weak  parts  of  such  a 
network  for  the  detection  and  identification  of  long -period  signals  from  explosions 
and  earthquakes.  The  evaluation  emphasizes  theoretical  characteristics  of  a 
network,  primarily  because  of  the  experimental  value  of  the  seismograph  systems 
and  a  consequent  limitation  of  joint  multiple  station  signal  observations.  Sufficient 
real  oata  are  given  to  demonstrate  the  theoretical  cases,  however,  and  later  data 
observations  can  be  expected  to  refine  the  estimc.ted  capabilities,  but  not  change 
the  characteristics  expressed  here  to  any  great  extent. 

Section  II  of  this  report  provides  the  theoretical  basis  for  estima¬ 
tion  of  signal  detection  capability  at  the  VLPE  stations  ?nd  the  observational  data 
supporting  the  estimate.  The  single  station  capability  is  then  merged  into  a  theo¬ 
retical  network  capability  model.  The  network  model,  which  may  include  contri¬ 
butions  from  the  long  period  arrays,  ALP  A,  NORSAR,  and  LASA,  is  then  used  to 
calculate  estimates  of  the  world- wide  surface  wave  detection  capability  of  networks 
of  stations  selected  from  the  current  and  projected  station  list.  Eoth  20  and  40 
second  Rayleigh  wave  detection  estimates  are  presented. 
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SECTION  II 


PRELIMINARY  ESTIMATES  OF  NETWORK  DETECTION  CAPABILITIES 
USING  ALPA,  NORSAR,  LASA,  AND  THE  VLPE  STATIONS 

A.  SINGLE  CHANNEL  DETECTION  MODEL 

1.  Basic  Model 

If  x  is  a  numerical  valued  random  phenomena  obeying  a  normal 
probability  law  with  parameters  (jz,  or)  representing  the  mean  (jz)  and  standard 
deviation  (<r)  of  the  distribution,  then  the  probability  that  the  random  value  x 
will  be  less  than  or  equal  to  X  is  given  by  (Parzen  1967): 

X  _  1  / ?  (  x 

r  1  1  r  1 '  a  ' 

P  I  x  <  X I  =  -  /  e  dx  (1) 

L  J  VTk*  J 

-CO 

We  will  apply  this  well-known  model  to  estimate  the  probability  of  long-period 
signal  detection  at  the  VLPE  stations.  Note  that  the  probability  that  a  randomly 
selected  value  will  be  less  than  or  equal  to  the  mean  of  the  distribution  is  given 
by  the  relationship  above  for  X  =  fi. ,  with  a  probability  of  0.  5,  and  the  probability 
that  it  will  be  greater  is  1.0  -  0.5  =  0.5,  since  the  function  is  continuous  { -co  +  <»}, 
The  point  of  interest  to  the  detection  problem  for  this  distribution  is  that,  given 
an  estimate  of  the  distribution  of  signal  and  noise  amplitudes,  a  standard  method 
is  available  for  estimating  the  probability  characteristics  of  detection.  Such  a 
method  may  be  applied  if  seismic  parameters  can  be  described  which  satisfy  the 
basic  model  requirements. 

2,  Seismic  Parameters  for  the  Defection  Model 

The  probability  field  for  signal  detection,  in  simple  form,  describes 
the  probability  that  the  signal  amplitude  (Ag)  from  a  seismic  event  wil1  exceed  the 
noise  amplitude  (A  )  so  that  a  "useful"  or  unbiased  measurement  of  the  signal  may 
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be  obtained  in  order  to  estimate  the  source  magnitude.  Both  signal  and  noise 
amplitudes  are  known  to  vary,  and  it  is  assumed  that  both  are  lognormal 
variables.  Wie  also  assume  that  the  signal  and  noise  are  independent  variables. 

Given  two  independent  variables  X ^  and  with  normal  dis¬ 
tributions  N(/ij,  <r  j)  and  respectively,  a  random  variable  describing 

the  distribution  of  the  difference  between  the  two  may  be  defined  which  suits  the 
requirements  for  determining  the  probability  that  a  random  signal  amplitude 

will  exceed  a  random  noise  amplitude.  Letting  X  be  the  distribution  of  the 

Z  1 

logarithms  of  signal  amplitude  N(  u  ,  a  )  and  X_  the  distribution  of  the  logar- 

^  S  S  L 

ithms  of  noise  amplitude  N(u  ,  a  )-  define  a  function  Y  : 

n  n 

V  =  X2-Xl  (2) 

According  to  distribution  theory  (Hogg  and  Craig,  1968),  the  new  random  vari- 

2  2 

able  Y  is  also  normally  distributed  N (u  -u  ,  a  +<r  ).  The  probability  that 

n  s  n  s 

Y  is  less  than  or  equal  to  zero  (i.  e. ,  that  the  signal  amplitude  is  greater  than 
or  equal  to  the  no* se  amplitude  calculated  from  the  basic  model  in  equation(l) 


Y  =  X2  -  Xj  (4) 


3.  Estimation  of  u  and  a 

s  s 

We  commonly  observe  that  the  estimate  of  surface  wave  magni¬ 
tude  (M  )  made  at  several  stations  shows  considerable  diffe*  ;nce  because  of 
s 
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radiation  pattern  at  the  source  or  path  dependent  signal  attenuation  effects. 

The  reported  magnitude  is  usually  an  average  of  the  observations,  and  an 

estimate  of  the  variance  of  the  M  observation  may  be  calculated.  The  M 

s  s 

value  here  is,  therefore,  assumed  here  to  be  a  normally  distributed  random 
variable  for  our  purposes,  with  a  mean  value  of  and  a  standard  deviation 
of  .  We  may  then  calculate  the  value  of  the  logarithm  t*f  the  signal  ampli¬ 
tude,  log  A^,  given  the  magnitude  of  the  event  and  the  distance  to  the  observing 
station. 

Assuming  a  simple  magnitude  model  for  surface  wave  observa¬ 
tions  near  20  seconds: 


where 


M  =  log  A/T  +  1. 12  +  logA  (+  «r  )  (5) 

S  —  M 

M  =  surface  wave  magnitude 

A  =  peak-to-peak  ground  displacement  in  millimicrons 
T  =  wave  period  in  seconds  for  A  measurement 
A  =  epicentral  distance  between  source  and  recording  station 
in  degrees  of  central  arc 

0 ^  =  standard  deviation  of  magnitude  data  based  upon 

estimated  from  several  stations  for  the  same  source  event 


If  A  is  expressed  in  kilometers,  T  =  20  seconds,  and  the  expected  value  of  A 
equal  to  Ag,  the  magnitude  becomes: 


M  =  log  A  -  2.  23  +  log  A,  (6) 

s  °  s  kms 

Then,  for  an  event  with  a  given  M  value,  the  amplitude  of  signal  expected  as 
some  distance  A  (kilometers)  is: 


log  A  =  M  +  2.  23  -  log  A  (7) 

s  s 

Assuming  that  all  variations  in  the  magnitude  estimate  are  in  the  log  A  term 
(i.  e.  ,  T  is  a  constant  20  seconds  and  distance  is  a  constant  for  a  given  observa¬ 
tion),  then  variation  in  the  logarithmic  term  represents  radiation  and  path  effects. 
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The  standard  deviation  of  the  logarithms  of  amplitudes  therefore  equals  the 
magnitude  standard  deviation,  or 

s  M  ' 

/*s  (±  *.)  =  Mg  *  2. 23  -  log  A  (8> 

4.  Estimation  of  fi  and  c 

n  n . 

Estimates  of  the  seismic  noise  amplitude  as  viewed  through  the 
seismograph  pass  band  may  be  calculated  from  direct  measurement  of  ampli¬ 
tude  on  visual  recordings  providing  proper  corrections  for  system  amplifica¬ 
tion  have  been  taken  in  order  to  determine  true  ground  motion.  For  this  pre¬ 
liminary  estimate  of  VLPE  noise,  the  Power  Spectral  Density  (PSD)  measure¬ 
ment  of  seismic  noise  has  been  obtained  from  a  number  of  selected  noise  sam¬ 
ples  of  about  one -hour  length  (actually  3584  seconds  length)  in  time  windows 
free  from  known  seismic  signals  as  an  alternate  method.  The  RMS  noise  (over 
a  given  frequency  band)  is  obtained  from  the  PSD  and  the  peak-to-peak  ampli¬ 
tude  is  estimated  by  multiplying  the  RMS  values  by  some  constant  for  an  esti¬ 
mate,  of  A^,  the  amplitude  of  noise. 

VLPE  digital  recording  systems  have  been  calibrated  to  provide 
a  measure  of  the  number  of  computer  counts  (C)  output  from  the  system  for  a 
given  input  of  ground  motion  in  millimicrons  (  I  ),  at  several  frequencies  within 
the  passband  of  interest.  The  PSD  provides  a  smoothed  estimate  of  the  density 
in  computer  counts  (C  /Af)  over  a  narrow  frequency  band  (Af)  centered  at  nAf 
steps  (n  =  0,  1,  2,  . ,  . ,  m),  with  the  final  value  m  equal  to  the  Nyquist  fre¬ 
quency  of  the  data  sample.  For  the  noise  estimate  here,  the  bands  ranging  from 
about  20  to  30  seconds  and  20  to  40  seconds  are  used.  The  RMS  noise  level  is 
simply  the  square  root  of  the  sum  of  the  PSD  (corrected  for  frequency  response) 
over  the  frequency  band  of  i.herest  (njAf  to  n^Af): 


sfioBse  amplitcsdes  ini  6lsc  2 0  to  30  second  period  range  are  pre- 
stmcffl  here  to  represent!  the  amplftode  of  noise  wbic&  sill  interfere  with  visual 
signal  detection  at  about  20  sec© nd  wave  period  in  se  j;rnagranas  nnproce.»sed 
except  for  the  seismograph  pass  ban  is  The  broader  band  of  2C-40  seconds  is 
used. to  estimate  interfering  noise  for  waves  of  about  40  second  period  since 
the  relatively  few  cycles  of  about  40  second  period  in  a  typically  dispersed  signal 
wave  train  may  be  more  difficult  to  discern  under  both  higher  and  lower  fre¬ 
quency  noise  conditions.  Arguments  can  certainly  be  made  thni  other  passbands 
would  be  more  representative,  but  these  bands  are  representative  of  those  which 
interfere  with  visual  analysis  procedures,  and  the  levels  are  representative  of 

t 

those  which  are  expected  to  interfere  with  signal  detection. 

Estimates  of  seismic  noise  characteristics  were  obtained  from 
seven  VLPE  sites  for  the  vertical  instrument  channel  by  converting  the  PSD 
measurements  to  RMS  ground  motion  (corrected  for  system  frequency  response) 
and  averaging  the  result  for  the  20  -  30  second  and  20  -  40  second  period  bands. 
Distribution  of  the  samples,  which  are  approximately  one  hour  in  length  and  free 
from  known  signals,  is  shown  in  Figure  II- 1  as  the  number  of  samples  taken  by 
Julian  day  (1971).  Table  II- 1  gives  the  average  of  observed  values  at  each  sta¬ 
tion,  with  data  for  the  Australian  location  estimated  on  the  basis  of  only  four 
observations  because  of  digital  recording  difficulties. 

5.  Interfering  Noise  and  Signal-to-Noise  Ratio 

If  the  RMS  amplitude  is  considered  as  a  measure  of  the  square 
root  of  the  mean  squared  amplitude  of  an  equivalent  steady  state  sinusoidal  sig¬ 
nal  at  the  frequency  of  interest  with  amplitude  A,  frequency  \Y,  and  initial  phase 
<p ,  then  instantaneous  amplitude  at  time  t  : 
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FIGURE  II- 1 

ONE-HOUR  NOISE  SAMPLING  AT  VLPE  SITES  FOR  NOISE  LEVEL  ESTIMATE 


TA3JLE  H-l 


MEAN  SE3SMBC  NOSKE  1  FOR  20-5®  SECOND  ME)  20-4®  SECOND 

BANDS,  AND  STAND aUdIj SVIATBON  OF  NOISE  IN  THESE  3*MB 

F«H  SEVEN  VLPE  SITES 


Station 

Noise 

m^*RMS 

20-50  Seconds 

•a 

(equivalent  in 
naagnitode) 
20-30  Seconds 

Australia 

(12-  1) 

(-23) 

Thailand 

12.6 

.19 

Alaska 

13.9 

.27 

Spain 

7.6 

.46 

Israel 

3.8 

.31 

Norway 

6.7 

.  17 

New  Jersey 

19-6 

.27 

mMRMS 

20-4C’  Seconds 

*h  (equivalent  in 
magnitude) 
20-40  Seconds 

Australia 

(15.0) 

(-28) 

Thailand 

12.8 

.22 

Alaska 

18.  1 

.29 

Spain 

8.0 

.42 

Israel 

4.  3 

.  33 

Norway 

7.4 

.  16 

New  Jersey 

27.6 

.25 

Given  a  signal  with  amplitude  A^=  A^fl_  x  2.  S28  at  'he  same 
frequency,  it  should  be  possible  to  slosr  a  signal  was  present  at  least  one -ball 
of  the  time  if  this  amplitude  equaled  the  noise  amplitude  {simply  dee  to  changes 
in  amplitude  and  phasing  caused  by  interference  of  the  two).  This  level  of  sigraal- 
to- noise  does  not  give  an  adequate  opportunity  for  measurement  of  an  unbiased 
signal,  so  some  greater  ratio  of  signal-to-noise  is  necessary.  An  arbitrary 
figure  of  something  greater  than  about  3.  0  x  ncise^^.  should  improve  the 
measurement  accuracy  significantly,  especially  considering  that  steady  state 
signals  and  noise  are  not  the  ordinary  conditions  for  detection. 

Measurement  of  the  maximum  peak -to -peak  noise  amplitude  and 
RMS  noise  estimate  at  the  same  wave  period  for  1200  second  long-period  noise 
samples  given  by  Harley  (1971)  indicated  that  a  factor  of  5  -  6  imes  noise 
represented  the  maximum  peak-to-peak  measurement  in  the  noise  sample.  A 
factor  of  4.  5  (13  dB)  was  used  in  that  study  to  estimate  improvement  in  detection 
of  chirp  filtering.  A  similar  study  by  Laster  (1970)  showed  that  Httle  or  no  in¬ 
fluence  of  noise  was  noted  in  chirp  filtered  data  when  used  to  estimate  magnitude 
if  the  y.t  ro-to-peak  maximum  signal  amplitudes  in  the  20  second  period  range 
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were  2-4  t fmes  the  noise, 


oiKse, 


BJMS 


esttrasate..  Even  when  zero-to-ptaS:  slgnal-to- 
rstios  were  I  -  2,  oasise  cttueaannanattiats  is  the  s%nci  resulted  Ira 


changes  of  else  magnitude  estimate  ©a  the  order  ©f  only  5.  1  emits. 

The  Implications  In  these  resells  lead  tt©  the  conrlrjslorr  Shat  z 

£ ac£©r  off  Aici  6  x  o©Ise_- for  detecting  a  “useful"*  rteaSs-tt©-cwea3r  signal  Is 

RM5  ~  ~ 

z  reasonable  estimate  ©ff  the  factor  for  use  In  the  theoretical  model  ©f  single 
station  detection.  Proof  ©f  the  estimate  acst  wait  for  additional  observational 
data,  which  Itself  will  provide  z  befe  >  .-r  estimate  of  the  long-term  influences  ©f 
EBon-selsmic  factors  as  well  as  regional  and  local  seismic  characteristics  of 
importance  t©  the  estimate. 


This  factor  of  6  Is  assumed  here  to  represent  both  a  conversion 
of  the  RAIS  noise  amplitude  to  "’Interfering  noise  amplitude”  for  signal  detection 
as  well  as  a  ""signal  to  noise  ratio®"  factor.  Considerable  argument  about  the 
proper  factor  for  use  has  been  presented  in  discussions  by  others  concerned 
with  the  purpose  and  use  of  th^~  results.  The  factor  should  be  doubled  according 
to  #feese  arguments,  i.  e. ,  a  factor  of  12  x  RMS  for  the  same  level  of  confidence 
in  detection.  This  appears  somewhat  overly  conservative  in  character  accord¬ 
ing  to  the  paragraph  above.  The  net  result  of  using  the  12  x  RliS  estimate  is  tc 
elevate  the  estimates  of  capability  presented  here  by  0.  3  M  units.  It  is  impor¬ 
tant  to  remember  that  the  6  x  RMS  (or  12  x  RMS)  represents  the  .  50  probable 
detectable  peak-to-peak  signal  amplitude  In  the  method  of  estimating  detection 
capability,  and  higher  or  lower  signal  amplitu  ies  will  have  greater  or  lesser 
probability  of  detection,  respectively. 

6.  Summary- 

Seismic  parameters  for  estimation  of  the  single  station  detection 
capability  using  the  normal  distribution  function  are: 
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+  o 
—  s 

=  M 

s 

+  2.23  -log  A, 

(14) 

^n 

±  a 

n 

=  log 
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RMS  n 

(15) 

siragle  sus&iaa  fflettecttsors,  casing  £  log  normal  dislirs&cairon: 
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NETWORK  DETECTION  MODEL 


I.  Basic  Model 


The  fundamental  network  model  used  in  this  study  was  developed 
under  contracts  AF3 3657- 12447  and  F33657-C-Q94I,  and  has  been  modified  for 
the  special  case  of  VLPE  data  where  the  method  of  estimating  the  noise  para¬ 
meter  is  consistent  with  the  signal  estimation  process. 

If  we  state  that  Aj,  ...»  are  the  observations  in  an  M- 
station  network,  with  probabilities  that  signal  will  exceed  noise  by  a  factor  of 
x  noted  as  P(Aj),  P(A^)t  . ..,  P(A^)  respectively,  then  the  probability  that 
some  number  of  stations  m  ,  P(m),  will  detect  may  be  calculated.  Following 
Parzen  (1967): 


so  ’  1 


XVI 

S1  =  £  P(\)  =  PfAj)  +  P(AZ)  FfAj^) 
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M  M 


s2=e  e  P(A  A  : )  =  PCA^)  *  PiAj 

kl  *2 

=  PtAjJP CA2)  4  PCAjJPt^)  4... 

(L  c,  ,  ail  non  repeating  2  station  products) 


M  Id 


S  = 
r 


*E  E  -E  piw . \> 


ki  k2 


(i.  e.  ,  all  non  repeating  r  station  products) 


SM  =  T>^AiA2» 


(Sr  is  the  summation  of  probabilities  over  all  possible  combinations, 

[  (r)  =  r!  (M-r) !  ’  r  =  2»  *  •  •  ’  M)j 


The  probability  that  exactly  m  stations  will  detect,  P(m)  , 


is  defined  by: 
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p(m, .  £  (;)  sr 


=  s 


(m+l\  /m+2\ 

m  j  Sm+1  +  (  m  ) 


t  *  •  «  » 
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and  the  probability  that  at  least  m  will  detect,  P{  >  m),  is 


(21) 


-  — ?.v  • 


VV.  v 


“W 


ji 


i 

i 


! 


I 


P«>m)  =  P(m)  -S-  P(m€i)  •£■  Ptn&SZ)  ■£•...  -6-  P(m4rs  =  M) 
M 

-  £  m  fc :  o 

r=m 


(22) 

(23) 


sra 


2.  Derived  and  Estimated  Input  Pari  meters 

Equivalent  ground  motion  for  PSD  (fc  x  1_  .  „)  and  <r 
2  2  1/2  KrdS  s 

(  #s  +  is  given  for  the  VLPE  stations  in  Table  II- 2,  with  6  x 

equaling  the  0.  50  detection  probability.  An  estimate  for  two  additional  VLPE 

sites  is  also  provided,  as  is  an  estimate  for  array  processed  capability  (time 

phased  and  sum)  for  ALP  A,  NORSAR,  and  LAS  A,  where  approximately  \/n 

improvement  over  single  instrument  is  assumed. 

3.  Theoretical  Network  Capability 


If  it  is  required  that  a  signal  amplitude  from  some  event  location 
exceed  the  noise  observed  (or  estimated)  at  the  stations,  and  fur  the:  that  this 
will  occur  at  some  minimum  number  of  stations,  a  probability  of  detection  by 
the  network  may  be  calculated.  Given  the  estimated  noise,  shown  in  Table  11-2, 
and  a  signal  amplitude  decay  with  distance  (equation  14),  then  the  individual 
station  probabilities  that  signal  exceeds  noise  are  calculated  (equation  17).  Thess< 
individual  probabilities  for  some  given  event  magnitude  and  location  are  then 
calculated  (equation  18)  and  the  network  probability  of  signal  detection  by  at 
least  m  of  M  stations  is  computed  (equation  22). 


4.  20-Second  Rayleigh  Waves 

The  calculations  above  were  completed  for  several  station  net¬ 
works  made  up  of  VLPE  and  Long  Period  Array  installations  using  a  hypothe¬ 
tical  grid  of  epicenters  as  sources.  The  results  indicate  that  it  is  highly  likely 
(.  90  probable  or  greater)  that  twc  or  more  stations  in  a  network  comprised  of 
the  seven  VLPE  stations  listed  in  the  upper  part  of  Table  II-2  will  detect  a  20- 
second  Rayleigh  wave  in  most  cf  Eurasia  if  the  source  magnitude  is  equivalent 


H-12 


TABLE  II -2 


STATION  LOCATIONS,  ESTIMATED  AMPLITUDE  (m#i,  p-p)  OF 
INTERFERING  SEISMIC  NOISE  (EQUALING  6  i:  noise  )  AND 
COMBINED  STANDARD  DEVIATION  OF  THE  LOGARITHMS 
OF  SIGNAL  AND  NOISE  AMPLITUDE 


to  at  least  an  3.  3  shallow  focus  continental  earthquake.  In  the  Kamchatka  - 

Kurile  region  (eastern  USSR),  an  M^=  2.4  source  magnitude  is  required,  and  in 

western  Eurasia  the  lowest  magnitude  detectable  under  these  criteria  is  about 

-  3.0.  Capability  for  tv/o  station  detection  is  about  M^=  3.2  in  central  Russia 

for  this  network,  and  aoout  M  =  3.  2  in  most  of  China. 

s 

If  the  VJLPE  stations  in  Norway  and  Alaska  are  deleted  from  the 
network  and  the  Long  Period  Arrays  in  Norway,  Montana,  and  Alaska  are  in¬ 
cluded,  20-second  Rayleigh  wave  detection  in  Eurasia  is  improved  by  about  0.  2 
magnitude  units  overall.  The  Kamchatka-Kurile  region  shows  slightly  greater 
improvement  to  a  two-station  capability  of  M  =  3.  1.  Inclusion  of  relatively  low 
quality  stations  in  Hawaii  and  Japan  reduce  the  detectable  magnitudes  0,  1 
units  in  the  Kamchatka-Kurile  region  and  in  eastern  China,  resulting  in  "thres¬ 
hold"  magnitudes  of  about  3.  0  for  highly  probable  two-station  Rayleigh  %vave 
detection. 

Conversion  of  to  equivalents  in  body  wave  magnitude  (m^)  of 
these  values  is  difficult  since  some  questions  have  been  presented  about  the 
relationship  between  the  two  at  low  magnitudes.  If  the  well  known  Gutenberg 
relationship  Mg=  1.  59  m^-  3.  97  is  used,  capability  for  the  seven-station  net¬ 
work  is  about  m,  =  4.  5  -  4.  6,  and  for  the  10  station  network  of  about  m ,  =  4.  3  - 
b  b 

4.4.  The  Gutenberg  relationship  is  usually  accepted  to  be  valid  above  rn^= 

4.  5  -  5.  0,  so  these  figures  must  be  considered  only  tentative. 

5.  40-Second  Rayleigh  Waves 

Examination  of  the  VLPE  data  indicates  that  surface  wave  am¬ 
plitude  attenuation  for  a  given  path  is  not  the  same  for  all  wave  periods.  In 
terms  of  signal  detectability,  it  is  especially  important  to  note  that  the  attenua¬ 
tion  of  40  second  Rayleigh  waves  is  greater  than  that  of  20  second  waves  on  the 
whole,  and  the  effectiveness  of  a  noise  "window"  at  the  longer  periods  is  there¬ 
fore  influenced.  Estimates  of  source  magnitude  are  also  affected,  but  insuf¬ 
ficient  information  has  been  published  to  estimate  this  factor  very  well. 
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Tryggvason  (1965)  showed  data  indicating  the  effective  ampli¬ 
tude  attenuation  (vertical  component  of  Rayleigh  motion)  as  a  function  of  wave 
period  and  distance  from  14  WWSSN  stations  recording  a  nuclear  explosion  on 
Novaya  Zemlya  in  1962.  The  attenuation  factor  given  by  Tryggvason  is  in  terms 
of  Q where  Q  is  the  effective  intrinsic  attenuation  factor  for  seismic  waves 
Tsai  and  Aki  (1970)  show  additional  data  in  the  20-  to  40-second  period  range 
which  essentially  agreed  with  Tryggvason's  result. 

The  point  of  interest  in  both  studies  for  application  here  is  that 
the  attenuation  of  40  second  waves  is,  for  average  world  paths,  apparently 
greater  than  that  for  20  second  waves.  The  interpretation  of  this  situation  is 
that  40  second  waves  penetrate  into  the  "low  velocity  zone"  of  the  earth's  mai.Te 
where  relatively  low  Q  exists,  and  40  second  amplitudes  are  reduced  at  a  great¬ 
er  rate  with  distance  than  the  20  second  amplitudes. 


If  we  assume  simple  geometrical  spreading  and  absorption  models 

for  the  20  and  40  second  waves,  an  estimate  of  the  relative  attenuation  influence 

on  magnitude  determinations  is  straightforward: 

7rfr 


where 


V 


is  the  amplitude  of  20  second  surface  waves 
r  kilometers  from  the  source 


(24) 


is  the  amplitude  of  20  second  surface  waves 
at  the  source 


rn  is  a  geometrical  spreading  in  kilometers 
f  is  surface  wave  frequency  (=  T  *  T=  wave  period) 


V  is  surface  wave  group  velocity 


Q20  is  the  effective  intrinsic  seismic  wave 
attenuation  factor  for  20  second  waves 
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and  for  40  second  waves  with  appropriate  subscript  changes: 


Ttix 

o  0  '°40v 

A40r  e 


(25) 


The  ratio  of  40  second  amplitudes  to  20  second  amplitudes  observed  at  a 
single  station,  letting  K  =  *s: 

qtv 


K20>r 


(26) 


assuming  that  the  amplitude  for  both  periods  is  the  same  at  the  source,  that 
the  geometrical  term  for  both  periods  is  the  same,  and  that  radiation  patterns 
are  independent  of  period. 


Figure  II-2  gives  ^ryggvason's  results  for  the  explosion  obser¬ 
vations  in  terms  of  K  vs.  wave  period.  Scaling  from  the  figure,  K  at  40  sec- 

—6  —  1  -6-1 
onds  is  about  160  x  10  km  ,  and  for  20  second  waves  about  80  x  10  km 

The  amplitude  ratio  computed  from  the  relationship  may  be  used  to  estimate 

the  amplitude  expected  at  40  seconds  given  a  20  second  measurement,  and  the 

magnitude  difference  expected  from  the  amplitude  difference  is  obtained  easily 

from  standard  magnitude  curves.  Equating  *o  estimate  through 

the  predicted  amplitude  difference  (r  =  A): 


M 


s40 


=  log  A  -  2.  53  +  log  A  +  .  35 


io"4a 


(27) 


The  difference  in  magnitudes  expected  between  20  and  40  second 
waves  due  to  difference  in  amplitude  attenuation  is  shown  as  a  function  of  dis¬ 
tance  for  six  of  the  VLPE  sites  in  Figures  II- 3(a)  and  II- 3(b).  Data  on  the  fig¬ 
ures  include  the  M  ..  and  M  _  „  calculated  from  observed  amplitudes  at  a  wave 
s40  s20 

period  of  about  20  and  40  seconds.  Sources  include  only  shallow  focus  earth¬ 
quakes  (33  kms  or  less  from  PDE  reports)  occuring  in  Eurasia  primarily,  and 


Presumed 

Explosion 


DISTANCE  (Degrees) 

FIGURE  II-  3  (b) 

DIFFERENCE  IN  MAGNITUDE  EXPECTED  DETWEEN  iiO  AND  40  SECOND  WAVES  DUE  TO 
DIFFERENCE  IN  AMPLITUDE  ATTENUATION 


several  known  or  presented  explosions  (frwlnsdtgg  one  in  Nevada  and  one  in 
the  Aleutian  Islands). 

While  the  amsant  of  data  available  is  Bruited,  these  resedas 
indicate  that  use  of  the  attenuation  difference  cooSd  be  a  valuable  tool  for  es¬ 
timating  source  magnitude  from  the  differ  eat  wave  periods.  Given  snore  acos- 
rate  hyyocenter  information,  some  of  the  scatter  in  magnitude  difference  can 
be  reduced  (deep  events  tend  to  result  in  about  equal  magnitude  estimates  from 
either  wave  period),  and  improved  estimates  of  signal  detectability  or  source 
influence  may  be  possible.  The  and  M^q  difference  reported  by  many 

authors  may  be  biased  by  the  distance  dependence  clearly  indicated  in  all  VjLPE 
data'sfeown  here,  and  some  correction  for  this  factor  must  be  taken  into  account 
to  evaluate  the  meaning  of  such  a  difference. 

In  addition  to  increasing  the  amount  of  data  for  purposes  of  im¬ 
proving  an  estimate  of  the  relative  attenuation  for  many  wave  paths,  some 
regionalization  of  effective  attenuation  will  no  doubt  be  necessary.  Paths  which 
show  about  the  same  relative  attenuation  for  several  events,  but  which  do  not 
necessarily  agree  with  the  predicted  relationship,  were  noted  within  the  data 
given  in  the  figures.  An  evaluation  of  these  regions  is  needed,  both  for  fuller 
understanding  of  the  attenuation  processes  and  for  delineati  ig  their  effect  upon 
detection  capability  and  scarce  mechanism  description. 

Using  the  same  grid  of  theoretica.  epicenters  as  was  used  for  the 

20-second  Rayleigh  wave  detection  estimate  results  in  a  slightly  low '  *■  detectic- 

threshold  at  40  seconds  when  the  epicenters  are  near  the  greatest  den  -  of  st; 

tion  positions.  A  rather  rapid  loso  of  this  extra  capability  with  distance  . 

because  of  the  additional  attenuation.  Most  of  eastern  Eurasia  (Europe)  is  ■  .  e* 

ed  by  highly  probable  (.  90),  2  or  more  station  detection  for  events  of  magnitudes 

Mg=  2.  8  -  2.  9  or  greater  using  the  basic  7-station  network.  The  addition  of  the 

3  large  arrays  depresses  the  threshold  here  by  0.  1  -  0,  2  M  units.  Along  the 

s 

USSR  -  China  border,  the  basic  network  has  calculated  capability  of  M  =  3.  1-3.2 
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Tbs  aa&fjgfina  of  arrays  again  adds  0.  5  -  0.2  M  crags  go  this  capaMligy,  Bov- 

s 

earing  tfee  g&resfeold  to  M^=  5.  0  -  3.  B. 

Capability  for  cSetectioa  of  40  ssctsad  Rayleigh  wares  for  sources 
in  the  Kamcfeatka-Kurile  region!  shows  the  effects  of  aflfterraaffga  strongly,  with 
tfee  basic  network  fearing  capability  of  M^=  3.  3  -  3.  4  and  addition  of  arrays 
dropping  tfee  threshold  fio  ahont  M  =  3.  0.  -Approximately  tfee  same  Jewels  are 
observed  in  central  and  soatfcerm  China.  Contribution  to  reducing  gfeis  ffert*$- 
feold  by  Hawaiian  and  Japanese  stations  is  only  a  local  effect  and  not  a  strong 
contribution  to  tfee  40-second  detention  in  a  network  sense. 

It  is  important  to  remember  tfeat  these  vaioes  are  for  surface 
wave  magnitudes  computed  at  40  seconds.  If  the  noise  levels  and  variances, 
signal  amplitudes,  and  decay  rates  {with  distance)  were  tfee  same  at  20  and  40 
seconds  then  these  values  would  teOJ  units  lower  than  tfee  equivalent  20- 
second  values.  Tfee  values  of  both  are  about  the  same  with  tfee  0.  3  units  offset 
in  the  40-second  estimate  because  higher  decay  rates  and  slightly  higher  noise 
levels  were  used.  The  estimates  suggest  that  detection  at  40  seconds  may  be 
slightly  lower  than  detection  at  20  seconds,  but  a  large  observational  data  base 
will  ba  required  to  make  a  quantitative  comparison. 
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NETWORK  PROCESSING  3¥  BEAMFOBMING  AND  ALSO  STACKING 

MATCHES  FILTER  OUTPUTS 


A-  GENERAL 

The  objectives  of  “time  delay  and  sum”  procedures  applied 
Co  the  VLPE  data  include  taking  advantage  of  the  opportunity  so  add  in-phase 
caprocessed  or  processed  signals  while  also  adding  randomly  phased  noise. 

Tfeea  r  etically  the  result  should  improve  she  signal  -So  -noise  ratio  by  a  factor 
of  VST  .  where  n  is  she  n- amber  of  inputs  So  the  summat2on,and  thereby  improve 
overall  capability  So  detect  signals  of  interest  for  further  analysis.  In  a 
practical  sense,  it  should  also  be  assumed  that  the  signal  amplitude  expected 
at  widely  spaced  single  sensors  will  have  come  reasonable  amplitude  in  com¬ 
parison  to  the  ncise. 

Long  period  signals  which  can  be  treated  in  this  way  include 
the  long -period  body  phases  P  and  S,  and  surface  wave  signals  processed  in 
some  way  so  that  in-phase  characteristics  could  be  obtained  for  useful 
summation  of  the  result.  Each  of  these  phases  has  certain  limitations  in 
both  physical  characteristics  and  possible  uses  for  discrimination  between 
earthquake  and  explosion  source  types.  Therefore  their  use  will  be  of 
interest  as  supporting  evidence  to  other  criteria  or  in  the  rare  situation 
where  these  signals  are  the  only  data  available  or  when  all  other  criteria  fail. 
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B.  JL0J8G  PEBJQD  BODY  PHASES 

Dteteotioa  of  Song-peri: sd  P  asd  S  has  several  factors  appro¬ 
priate  Co  the  discrimination  problem.  One  of  these  factors  is  that  the  dMf- 
v  recce  in  energy  expected  between  tire  long -period  phases  and  higher  fre- 
sjueecy  f  short -period)  body  waves  shocM  he  oa  the  order  of  Shat  indicated 
by  the  ‘"standard  discriminant05  of  la  other  words,  an  m^  or 

most  important  when  the  surface  ware  is  not  detectable.  As  another  example, 
while  demonstrati '‘as  of  significant  shear  energy  does  not  exclude  the  possibility 
of  tectonic  shear  energy  release  in  the  vicinity  of  an  explosion,  a  total  lack  of 
shear  energy  might  suggest  a  very  small  effective  source  volume  of  "explosion- 
like’"'  scarce  mechanism. 

Signal  detection  capability  for  short -period  P-waves  is  gen¬ 
erally  better  than  that  for  the  long -period  body  phases.  For  this  reason, 
hypocenter  parameters  are  available  which  provide  a  reference  for  reasonably 
accurate  estimation  of  relative  time  delays  for  the  long-period  body  phases 
within  a  world-wide  or  continental-dimensioned  network  of  recording  stations. 
Some  loss  in  the  theoretical  \/n~  improvement  can  be  expected  due  to  poss¬ 
ible  radiation  patterns  and  selective  travelpath  effects  (such  as  frequency 
dependent  energy  absorption),  but  if  some  library  of  information  for  regions 
of  interest  can  be  accumulated,  the  effect  of  such  factors  can  be  minimized. 

When  more  multiple  station  digital  data  become  available, 
multiple  observations  of  long -period  P  and  S  will  permit  an  evaluation  of 
the  beam-forming  potentials  for  the  VLPE  network  with  the  following  goals: 


:  m^  might  be  developed.  The  utility  of  such 


becomes 


•  DesnoESStraEe  Vm  Bnapjrs»Tenmcffltt  or  lack  of  improvemenit 

•  EstaMss^  SEEDglc  a»<ffwfl  multiple  sttattio!;  sigmal  detections 
StmiL.  for  She  long-pciriod  body  phases 

•  Demonstrate  possible  discriminants  where  the  standard 

M  discriminant  cannot  be  obtained. 
s  b 

•  Accumulate  long-period  body  wave  characteristics  for 
regionalization  of  source  and  path  effects. 

An  example  of  the  beam-forming  approach  is  shown  below 
for  a  three-station  beam.  The  polarity  radiation  problem  for  P  probably 
does  mot  enter  in  this  example,  but  future  plans  include  an  evaluation  of 
the  possibility  of  improving  P  detectability  and  minimizing  the  polarity 
effect  by  forming  a  Z  x  R  beam  prior  to  time  delay  and  summation.  Z  x  R 
denotes  a  point-by-point  multiplication  of  simultaneous  vertical  and  radial 
data  to  enhance  in-phase  components,  such  as  P  motion,  at  the  expense  of  the 
out  of  phase  vertical-horizonal  motions  in  a  plane  radial  to  the  source.  The 
approach  is  not  useful  for  S  since  the  particle  motion  of  S  becomes  very 
complex  after  the  first  observed  cycle  or  half-cycle  in  most  cases. 

In  this  preliminary  test  of  the  LXTRAN  beamforming  program 

using  VL.PE  data,  a  presumed  explosion  at  49.7  N,  78- 2E,  m  =  5.8  was 

b 

selected  from  the  limited  multiple  station  observations  available.  This 
selection  was  made  in  order  to  eliminate  the  problem  of  signal  polarity 
which  might  be  a  problem  with  an  earthquake  source.  Delay  times  were 
computed  from  J-B  travel  times  contained  within  the  program,  then  applied 
to  the  data  traces  from  Thailand  and  Alaska.  The  third  station.  New  Jersey, 
was  used  as  a  reference  trace.  All  other  data  are  delayed  to  match  the 
P  arrival  times  with  the  expected  P  at  New  Jersey,  and  the  summation  trace 
time  reference  is  the  same  as  the  New  Jersey  traces.  Results  are  shown  in 
Figure  III— 1- 


Seamforming  was  unsuccessful  in  extracting  the  P  »«ve,  which 
is  not  considered  too  surprising  in  view  of  the  total  absence  of  visual  P  motion 
on  any  of  the  traces  and  the  very  few  stations  available  for  the  beam.  The 
amplitude  of  P  from  the  suspected  explosion  is  very  low,  probably  less  than 
1/ 10  the  amplitude  of  Rayleigh  observed  at  the  Thailand  and  Alaska  stations, 
suggesting  that  the  P  summation  will  be  useful  for  discrimination  primarily 
in  an  extreme  case  of  interference  of  the  higher  amplitude  Rayleigh  motion. 


C.  TIME  DELAYS  AND  STACKED  MATCHED  FILTER  OUTPUTS 

The  dispersive  character  of  surface  waves  and  the  difference  in 
dispersion  for  different  travel  paths  make  direct  "delay  -  and  -  sum"  tech¬ 
niques  impractical  for  these  propagation  modes.  For  this  reason,  delay  and 
sum  would  have  to  be  performed  on  matched  filter  outputs  which  ideally  com¬ 
press  the  surface  wave  signal  into  a  single  pulse-like  signaL 

Two  major  matched  filter  techniques  are  being  investigated:  (1) 
the  Master  Event  method,  where  surface  wave  signals  from  high  signal-to- 
noise  ratio  recordings  are  used  as  a  filter  to  emphasize  low-level  signals  from 
the  same  source  region,  and  (2),  Chirp  Filters,  which  attempt  to  provide  an 
inverse  of  the  surface  wave  dispersion  characteristic  of  some  particular  travel 
path.  Both  techniques  should  provide  a  strong  filter  output  when  the  filter 
"wavexorm"  matches  the  signal  waveform,  and  this  output  may  be  correct 
onougb  to  attempt  the  delay  and  sum  technique.  Close  match  of  signal  and  filter 
characteristics  is  necessary  to  obtain  the  compactness  needed.  A  large  library 
of  filters  is  necessary  because  of  the  strong  differences  of  surface  wave  dis¬ 
persion  among  the  possible  world  travel  paths. 
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Figures  IH-2  and  III -3  demonstrate  the  difficulties  related  to 
beamsteer  in  the  case  of  Marched  Filter  stacking.  The  time  delayed  summation 
of  signals  recorded  at  ALPAaod  NORSAR  from  an  event  in  eastern  Russia  (53. 4N, 
120.  3E)  were  filtered  with  a  master  waveform  and  chirp  filters.  The  results  show 
that  the  compact  signal  needed  is  not  formed  very  satisfactorily  and  if  the  beams 
were  stacked  according  to  a  ?5  second  Rayleigh  (LR)  velocity  of  3.  5  km/ sec  or 
Love  wave  (1X1)  velocity  of  about  3.  8  km/  sec  for  the  same  period  (noted  by  an 
arrow  on  the  figures),  not  too  much  information  would  be  obtained.  The  develop¬ 
ment  of  more  accurate  chirps  or  selection  of  Master  waveforms  which  more 
closely  model  the  signals  are  obviously  needed  in  this  instance,  again  indicating 
the  need  for  a  library  of  waveforms. 

D.  INTERFERING  EVENTS 

The  interfering  (or  mixed)  event  problem  is  significant  on  long 
period  surveillance;  Table  HI-1  (Mack,  1971)  shows  the  percentage  of  events  in¬ 
terfered  with  at  LASA  over  a  50 -day  period.  Generally  the  percentage  of  mixed 
events  is  about  15%  at  a  single  station,  and  is  dependent  on  magnitude  (the  problem 
is  more  severe  for  smaller  events). 

The  availability  of  multi-station  data  should  reduce  the  severity  of 
the  problem  somewhat,  because  the  distribution  of  individual  stations  increases 
the  probability  that  the  desired  event  will  be  "clearly"  recorded  at  some  station 
in  the  network.  However,  we  wish  to  obtain  measureable  signals  from  as  many 
locations  as  possible  so  that  maximum  information  about  the  signal  source  is  avail¬ 
able.  Techniques  for  separating  these  mixed  signals  are  being  investigated  for 
specific  problem  conditions,  since  the  ability  to  separate  improves  the  overall 
network  aspects  of  detection  and  discrimination  capability. 

In  some  instances  it  is  possible  to  separate  mixed  events  at  arrays. 

It  has  been  established  that  if  the  two  events  have  on  the  order  of  30°  or  more  azi¬ 
muthal  separation,  array  processing  techniques  can  be  used  to  extract  the  desired 
event  from  the  overlapping  event.  The  best  capability  is  no  more  than  one  magni¬ 
tude  unit  (i.  e.  ,  an  M  =  4.  0  event  could  be  recovered  from  an  M  =  5.  0  event). 

s  s 


Predicted  Arrival  Time,  25  Second  Energy 
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ING  FILTER  OUTPUTS  (POOR) 


TABLE  III-l 

PERCENTAGE  OF  MIXED*  EVENTS  AT  LAS  A  (MACK,  1971) 


Magnitude 

^  3.  9 

4.0  -  4.4 

4. 5  -  4. 9 

5.  0  -  5.5 

>5.5 

LAS  A 

Number  of  Events 

357 

239 

105 

52 

27 

%  Mixed 

19.9 

17.2 

19.1 

11.6 

7.4 

*  A  mixed  event  is  the  recording  of  multiple  signal  arrivals  from  different 
sources  on  the  seismograms  which  cannot  be  interpreted  with  enough  resol¬ 
ution  to  provide  useful  information  from  the  recorded  data. 
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For  mixed  events  having  about  the  same  azimuth,  array  process¬ 
ing  cannot  be  used.  However,  the  complex  cepstrum  technique  (Linville,  1971) 
can  be  used  for  this  case.  This  technique  also  could  be  applied  to  a  single  station 
when  the  mixed  events  have  different  azimuths. 

The  capability  of  the  complex  cepstrum  technique  to  separate  mixed 
events  has  not  been  fully  tested.  Figure  III-4  shows  results  for  a  test  case  for  two 
events  from  essentially  the  same  epicenter  (on  the  mid-Atlantic  ridge)  recorded  at 
Ogdensburg.  The  first  event  is  about  one-half  the  amplitude  of  the  second  (over¬ 
lapping)  event;  the  complex  cepstrum  separates  these  two  events  very  well. 

To  test  the  resolution  of  the  technique,  the  estimate  of  first  event 
was  scaled  down  by  a  factor  of  two  and  four,  and  added  back  to  the  second  event 
with  the  original  time  separation  (108  seconds).  Note  that  the  first  event  now  is 
one-fourth  (0.6  Mg  units)  and  one-eighth  (0.  9  Mg  units)  the  size  of  the  second 
(overlapping)  event.  Figures  III-5  and  III-6  show  that  separation  is  achieved  by 
applying  the  complex  cepstrum  technique  in  both  cases.  Note  that  when  the  first 
event  is  one-eighth  the  size  of  the  second  event,  it  is  difficult  to  see  it  in  the  mix¬ 
ed  waveform,  and  impossible  to  get  an  accurate  M  estimate. 

s 

In  this  case  the  complex  cepstrum  technique  recovered  the  desired 
event  when  it  was  0.  9  Mg  units  smaller  than  the  overlapping  event;  this  is  essen¬ 
tially  as  good  as  the  array  processing  capability.  It  should  be  noted  that  the  ex¬ 
ample  used,  while  constituting  real  data,  probably  represents  close  to  the  best 
case  for  this  technique  because  the  two  events  were  simple  (i.  e. ,  no  multi- 
pathing)  and  had  essentially  the  same  spectrum.  For  more  complicated  cases  the 
performance  probability  would  not  be  as  good,  but  still  should  offer  some  capabil¬ 
ity  to  separate  mixed  events.  The  technique  has  excellent  potential  because: 

•  I*  provided  separation  capability  for  mixed  events  having  similar 

epicenters  and,  in  fact,  probably  works  best  for  this  situation. 

Note  that  array  processing  is  ineffective  here. 

•  It  can  be  applied  at  a  single  station. 
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FIGURE  III-4 


RECOVERY  OF  A  LOW  LEVEL  SIGNAL  FROM  A  MIXED  EVENT  WITH 
SIGNAL-  TO-SIGNAL  RATIO’OF  ABOUT  2:1  (MAXIMUM  AMPLITUDE) 
BY  THE  COMPLEX  CEPSTRUM  TECHNIQUE  APPLIED  TO  SINGLE 
STATION  DATA.  MID- ATLANTIC  RIDGE  EVENTS  RECORDED  ON 
THE  OGDENSBURG  VLPE  SEISMOGRAPH  . 
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FIGURE  III -5 

RECOVERY  OF  A  LOW  LEVEL  SIGNAL  FROM  A  MIXED  EVENT  WITH 
SIGNAL -TO -SIGNAL  RATIO  OF  ABOUT  4:1  (MAXIMUM  AMPLITUDE) 
BY  THE  COMPLEX  CEPSTRUM  TECHNIQUE  APPLIED  TO  SINGLE 
STATION  DATA.  MID-ATLANTIC  RIDGE  EVENTS  RECORDED  ON 
THE  OGDENSBURG  VLPE  SEISMOGRAPH 
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FIGURE  III- 6 


RECOVERY  OF  A  LOW  LEVEL  SIGNAL  FROM  A  MIXED  EVENT  WITH 
SIGNAL -TO -SIGNAL  RATIO  OF  ABOUT  8:1  (MAXIMUM  AMPLITUDE) 
BY  THE  COMPLEX  CEPSTRUM  TECHNIQUE  APPLIED  TO  SINGLE 
STATION  DATA.  MID-ATLANTIC  RIDGE  EVENTS  RECORDED  ON 
THE  OGDENSBURG  VLPE  SEISMOGRAPH  ' 
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Mai ?i- station  digital  data  have  not  been  available  in  sufficient 
quantity  to  allow  os  to  quantize  tie  redaction  in  tie  somber  of  mined  events  pro¬ 
vided  by  tie  VJLPE  network.  We  plan  to  investigate  tins  in  tie  future  by  following 
steps  outlined  below: 

•  Select  a  continuous  time  period  of  one  to  two  mentis  in  duration 
(a  period  when  the  ViPE  network  has  many  stations  operational). 

•  Obtain  a  cc-mpr  Aensive  event  list  (Note  that  tie  list  currently  be¬ 
ing  generated  by  Lincoln  Laboratories  for  the  February-March 
time  period  would  be  suitable). 

•  For  each  event  and  each  station,  tabulate  the  percentage  of  mixed 
events.  In  our  application  we  would  limit  the  desired  events  to 
those  ficm  the  Sino -Soviet  area. 

•  Combine  the  station  data  to  produce  additional  listings  of  percent¬ 
age  mixed  at  all  stations,  all  but  one  station,  and  all  but  two  sta¬ 
tions. 

•  Compare  ‘.the  single  station  and  network  "mixing  levels"  as  a  func¬ 
tion  of  magnitude. 

The  result  of  this  study  will  be  to  provide  a  quantitative  estimate  cf  hew  much  the 
network  reduces  the  mixed  event  problem,  the  severity  of  the  problem  remaining, 
the  critical  stations  in  the  network  (these  should  be  determined  by  their  related 
location  with  respect  to  areas  of  high  seismicity  and  the  Sino-Soviet  area),  and 
the  improvement  achieved  (if  any)  by  increasing  the  number  of  stations. 

Included  in  the  tabulation  of  mixed  events  described  earlier  will  be 
a  measurement  of  the  amplitude  of  the  interfering  event  during  the  expected  arrival 
time  of  the  desired  event.  The  expected  amplitude  of  the  desired  event  will  be  ob¬ 
tained  (either  from  Ms  measurements  at  stations  where  the  desired  event  is  cl^-r- 

ly  recorded  or  by  converting  event  m  values  to  M  )  and  the  amplitude  ratio 

b  s 

A  A  will  be  calculated: 
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where; 


AA  = 


Aj  is  the  amplitude  of  the  interfering  event 

Ajj  is  the  estimate  of  the  amplitude  of  tie  desirec  event. 


Tie  quantity  A  A  provides  an  estimate  of  tie  *-eIative  size  of  tie  tiro  events.  7L, 
using  typical  values  for  array  processing  and/or  complex  cepstrxim  gains,  tie 
redaction  in  tie  percentage  of  mixed  events  which  could  be  achieved  will  be  es¬ 
timated.  The  residual  percentages  of  events  mixed  at  all,  all  bat  one,  and  all 
but  two  stations  provide  estimates  of  tie  size  of  tie  "anresolvable"  mixed  event 
problem. 


A  special  case  that  tie  above  analysis  may  cot  account  for  is  tie 
very  large  earthquake,  which  can  generate  large  amplitude  surface  waves  lasting 
several  hours.  Events  will  be  recorded  during  this  time  period  and 
many  will  not  be  recovered,  even  when  processing  gains  are  taken  into  account. 
The  mixed  event  problem  in  this  case  can  be  quantized  as  follows: 


•  Measure  coda  amplitude  decay  rates  for  several  large  earthquakes 
(Mg  -  6.  0)  at  several  stations  and  obtain  an  average  decay  rate. 
Both  network  and  station  averages  will  be  obtained. 

•  Determine  the  time  window  when  the  events  sought  will  be  unre¬ 
coverable  at  all  stations,  taking  processing  gains  into  account.  The 

time  window  will  be  magnitude  (M  )  dependent. 

s 

•  From  seismicity  data,  estimate  the  number  of  events  in  a  given 
magnitude  range  which  would  occur  during  the  time  period. 

•  -Again,  from  ~?,ismicity  data,  estimate  the  number  of  events  with 
Mg  -  6.  0  per  year. 

This  will  provide  an  estimate  of  the  number  of  unrecoverable  events  per  year  due 
to  very  large  earthquakes.  Depending  on  the  data  base  selected  in  the  analysis 
discussed  previously,  this  number  may  have  to  be  added  to  that  obtained  in  that 
study. 
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Another  special  case  is  tie  major  aftershock  sequence  associated 
with  a  large  earthquake  in  the  earthquake  swarm.  Both  sequences  typically  in¬ 
volve  several  days  time,  daring  «Mc!s  a  very  large  percentage  of  desired  events 
will  be  mixed.  Because  of  tne  complexity  of  these  sequences  a  detailed  study  of 
the  mixed  event  problem  would  be  difficult,  however,  the  following  simple  approach 
should  define  the  'worst  case "  situation. 

•  Determine  the  time  in  the  sequence  when  the  events  are  no  longer 
"continuous”  (i.  e.  ,  the  time  separation  between  events  exceeds 
event  duration  -  about  30  minutes  to  one  hour,  depending  on  the 
event  size). 

•  From  the  magnitudes  of  the  aftershocks  preceding  this  time,  es¬ 
timate  the  minimum  recoverable  magnitude, taking  processing  into 
account.  Note  that  this  magnitude  probably  would  be  time  depen¬ 
dent  and,  in  the  case  of  the  aftershock  sequence,  the  first  several 
hours  would  have  to  be  excluded  because  the  main  shock  would 
dominate. 

•  Again,  from  seisruicity  data,  estimate  the  number  of  mixed  events 
expected  during  the  time  period  and  the  number  of  swarms  and 
aftershock  sequences  per  year. 

This  procedure  would  provide  a  roug.i  estimate  of  the  number  of  mixed  events / year 
due  to  aftershock  sequences  in  swarms.  Again,  depending  on  the  original  data  base 
this  number  may  have  to  be  added  to  previous  values  to  obtain  a  complete  estimate 
of  the  mixed  event  problem. 

We  plan  to  look  at  these  special  cases  during  the  coming  year 

also. 
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SECTION  IV 
CONCLUSIONS 


The  VLPE  network  has  the  theoretical  capability  for  detection  of 
shallow  focus  continental  earthquakes  in  Eurasia  at  =  2.9  -  3.  0  when  in  opera¬ 
tional  status  and  in  combination  with  the  large  arrays.  The  theoretical  values  of 
the  estimate  requires  that  a  systematic  evaluation  of  actual  signal  detections  be 
reserved  as  the  final  demonstration  of  capability.  The  preliminary  empirical 
data  show  that  no  serious  error  in  estimation  is  likely,  but  some  refinement  will 
very  likely  be  possible.  This  is  particularly  true  in  the  case  of  40  second  surface 
waves,  which  are  attenuated  at  a  different  rate  than  the  20  second  signals. 

Because  of  limited  observational  data,  a  true  picture  of  die  net¬ 
work  contribution  to  mixed  event  and  radiation  pattern  problems  cannot  be  clearly 
demonstrated  at  this  time.  Beamforming  for  long  period  P  and  S  signal  detection 
may  have  utility  in  special  cases,  but,  these  are  not  likely  to  contribute  to  identi¬ 
fication  of  sources  as  routine  discriminants.  The  power  of  the  complex  cepstrum 
method  at  the  single  station,  when  fully  evaluated,  may  contribute  to  the  network 
aspect  very  significantly  in  terms  of  mixed  events.  The  approach  discussed  for 
description  of  a  set  of  mixed  event  parameters  in  Section  HI  can  provide  an  im¬ 
portant  measure  of  the  potential  for  this  approach.  Utility  of  matched  filter  stacks 
for  overall  reduction  of  the  detection  threshold  will  be  marginal  until  a  suitable 
library  of  filters  can  be  developed.  This  factor  influences  the  ability  of  a 
limited  number  of  stations  to  observe  the  radiation  patterns  with  enough  resolu¬ 
tion  to  describe  the  problem  in  any  detail. 
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